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ABSTRACT 


Gamma-ray bursts (GRBs) are classified as long and short events. Long GRBs (LGRBs) are associated with 
the end states of very massive stars, while short GRBs (SGRBs) are linked to the merger of compact objects. 
GRB 2008264 was a peculiar event, because by definition it was a SGRB, with a rest-frame duration of ~ 0.5 
s. However, this event was energetic and soft, which is consistent with LGRBs. The relatively low redshift 
(z = 0.7486) motivated a comprehensive, multi-wavelength follow-up campaign to characterize its host, search 
for a possible associated supernova (SN), and thus understand the origin of this burst. To this aim we obtained a 
combination of deep near-infrared (NIR) and optical imaging together with spectroscopy. Our analysis reveals 
an optical and NIR bump in the light curve whose luminosity and evolution is in agreement with several LGRB- 
SNe. Analysis of the prompt GRB shows that this event follows the Ep, — Eiso relation found for LGRBs. The 
host galaxy is a low-mass star-forming galaxy, typical for LGRBs, but with one of the highest star-formation 
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* Based on data obtained with the LBT programs: LBT-2019B-18 (PI A. 
Rossi), DDT-2019B-9 (PI P. D’ Avanzo), and LD-2020B-0100 (PI. B. Roth- 
berg); and the TNG programme A41. TACIS (PI V. D' Elia). 
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rates (SFR), especially with respect to its mass (log M./Mo = 8.6, SFR ~ 4.0 Mo/yr). We conclude that GRB 
200826A is a typical collapsar event in the low tail of the duration distribution of LGRBs. 

These findings support theoretical predictions that events produced by collapsars can be as short as 0.5 s in 
the host frame and further confirm that duration alone is not an efficient discriminator for the progenitor class 


of a GRB. 


Keywords: Gamma-ray bursts (629) — Supernovae (304) — Core-collapse supernovae (1668) 


1. INTRODUCTION 


Traditionally, gamma-ray bursts (GRBs) are divided in two 
classes based on their duration and their spectral hardness in 
the gamma energy range. The majority of observed GRBs 
have durations longer than 2 s, a soft high-energy spectrum 
and are termed “long GRBs” (LGRBs), as opposed to “short 
GRBs” (SGRBs) with a usually sub-second duration and a 
harder high-energy spectrum (Mazets et al. 1981; Kouve- 
liotou et al. 1993). However, there is a substantial over- 
lap between the distributions of the durations of the two 
GRB classes (Kouveliotou et al. 1993, Bromberg et al. 2013, 
Horvath & Toth 2016). 

Long GRBs are typically associated with highly energetic 
broad-lined type Ic supernovae (SNe) and have therefore 
been associated confidently with the deaths of very mas- 
sive stars (e.g., Woosley & Bloom 2006; Cano et al. 2017a). 
On the other hand, short GRBs are thought to be connected 
with the merger of compact objects, as was firmly demon- 
strated in spectacular fashion by the SGRB 170817A which 
was accompanied by the binary neutron star (NS) merger 
gravitational-wave source GW170817 detected by the LIGO 
and Virgo interferometers as well as by the kilonova (KN) 
AT2017gfo (e.g., Abbott et al. 2017a,b,c; Pian et al. 2017; 
Smartt et al. 2017). 

Improving our understanding of the connection between 
LGRBs and SNe has been one of the most intense areas of re- 
search in the field of GRBs for over 20 years. More than 1400 
bursts have been discovered by the Neil Gehrels Swift Obser- 
vatory (Swift hereafter) in the last 16 years, among which 
only 40 — 50 GRBs have an associated SN identified by the 
late bumps in their optical afterglow light curves, and to date, 
just 28 have been spectroscopically confirmed (Rossi et al. 
2022, in prep., Izzo et al. 2019; Melandri et al. 2021, Cano 
et al. 2017a and references therein, Klose et al. 2019; Cano 
et al. 2017b; Ashall et al. 2019a; Melandri et al. 2019; Hu 
et al. 2021). The burst duration for all these events lasted for 
more than 2 s, therefore they are considered LGRBs. While 
these results show a simple and clear correlation between 
SNe and LGRBs, the reality may be more murky. 

A first break in the LGRB-SN association was already 
found in 2006 with the discovery of the LGRBs 060505 


and 060614!. Their relatively low redshifts (z = 0.089 and 
z = 0.125, respectively) allowed for the search of the possi- 
bly associated SNe, but despite achieving deep photometric 
limits no SN component as luminous as those observed in 
the case of other long GRBs was detected (Della Valle et al. 
2006; Gal-Yam et al. 2006; Fynbo et al. 2006; Ofek et al. 
2007; McBreen et al. 2008; Xu et al. 2009)?. Later, LGRB 
111005A at z = 0.0133 (Tanga et al. 2018; Michałowski et al. 
2018) also joined this group. On the other hand, in the case 
of GRB 040924, classified as short for its rest-frame dura- 
tion (~ 1 s), an associated SN was detected (Wiersema et al. 
2008; Soderberg et al. 2006) thus identifying its origin as be- 
ing the collapse of a massive star. Another peculiar case is 
GRB 090426, which has a duration shorter than 2 s, but with 
a high-energy spectrum and afterglow brightness typical of 
long GRBs (Antonelli et al. 2009; Thóne et al. 2011; Nicuesa 
Guelbenzu et al. 2011, 2012; Kann et al. 2022, in prep.). In- 
deed, although classifying a burst by its duration as either 
long or short is still the most simple and adopted approach, 
the above results show that this is not enough to reveal the 
nature of the GRB progenitor. Even considering the spectral 
hardness often may not help in the classification (e.g., Mi- 
naev & Pozanenko 2020a; Agüí Fernández et al. 2021), and 
thus understanding the origin of the burst. Indeed, the most 
secure way is to determine whether the event is associated 
with a supernova (SN) and is a LGRB, or with a KN and is a 
SGRB. 

We here present follow-up observations and analysis of 
GRB 2008264, including the search for an associated SN 
and characterization of the host galaxy. This very intense, 
but short-duration burst was detected by several space-based 
GRB detectors (see Hurley et al. 2020), including Konus- 
Wind (Ridnaia et al. 2020; Svinkin et al. 2020), Fermi/GBM 
(Mangan et al. 2020), AGILE/MCAL (Pittori et al. 2020) and 
AstroSat/CZTI (Gupta et al. 2020). Stemming from its ob- 
served short duration (GBM Too = 1.14 s, Mangan et al. 


! GRBs 060505 and 060614 have durations (Too) of 4 + 1s (Hullinger et al. 
2006) and ~ 102 + 5 s (Barthelmy et al. 2006), respectively. 

? More recently, KNe have been claimed to have been detected in the after- 
glows of GRB 060614 (Yang et al. 2015; Jin et al. 2015) and GRB 060505 
(Jin et al. 2021). This would imply they are actually merger events (see 
Kann et al. 2011 for more discussion). 
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2020; Konus-Wind Too = 0.772 s, Svinkin et al. 2020), it 
was initially thought to be a genuine SGRB. The ZTF survey 
project (Bellm et al. 2019) launched an observing campaign 
to search for the optical counterpart (Sagues Carracedo et al. 
2020) within the large error circle of Fermi/GBM. Among the 
several candidates detected, ZTF20abwysqy was confirmed 
to be the afterglow of GRB 2008264 (Ahumada et al. 2020a; 
Belkin et al. 2020a; Dichiara et al. 2020). An X-ray after- 
glow was also detected by Swift/XRT (D' Ai et al. 2020). The 
GRB spectroscopic redshift (z = 0.748) was obtained thanks 
to LBT/MODS observations (Rothberg et al. 2020a). The 
provided distance was used to calculate the real luminosities 
and energy values, which showed that GRB 200826A was a 
very energetic event with an isotropic energy release and a 
peak energy more consistent with classical LGRBs than with 
typical SGRBs (Svinkin et al. 2020). A SN bump in the af- 
terglow light curve was finally reported by Ahumada et al. 
(2020c). Additional optical observations of the afterglow and 
the SN have been reported in Ahumada et al. (2021); Zhang 
et al. (2021). 

Throughout this work, the flux density of the afterglow is 
described as F,(t) œ v£, and a ACDM world model with 
Quy = 0.308, O4 = 0.692, and Ho = 67.8 km s^! Mpc^! 
(Planck Collaboration et al. 2016) has been assumed. All 
data are in observer frame, unless differently specified. 


2. OBSERVATIONS AND DATA ANALYSIS 
2.1. Optical Spectroscopy 


Optical spectroscopy of GRB 200826A was obtained 
with the Large Binocular Telescope (LBT) on Mt. Gra- 
ham, Arizona, USA, on August 28, 2020 at UT 09:06:45 
(MJD 59089.379692) using both optical Multi-Object Dou- 
ble Spectrographs (MODS, Pogge et al. 2010). Each MODS 
was configured using a 1" wide slit (R ~ 1500). MODS-1 was 
configured with the red grating (0.58 um — 1.0 um coverage) 
due to issues at the time with the blue channel, and MODS-2 
was configured with the dual grating (0.32 um — 1.05 um cov- 
erage). Observations consisted of 4 x 900 second exposures 
for the three channels, which yielded an effective exposure 
time of 7200 seconds for 0.58 um — 1.0 um (the red chan- 
nels). A position angle (PA) = 17.0? was used in order to 
cover both the reported location of the GRB 2008264 after- 
glow and a nearby galaxy 15" NE. Observations were ob- 
tained at or near transit with an effective airmass of 1.0005. 
The seeing for the spectroscopic observations ranged from 
077 — 1*0 and conditions were photometric. 

The MODS data were reduced first with modsCCDRed 
version 2.04 to remove the bias and flat-field the data. Next, 
custom IRAF scripts were used to: extract along the cen- 
tral slit using a stellar trace; wavelength calibrate and rectify 
the tilt in both X and Y directions using arc-lamp lines; flux 
calibrate the data using the spectro-photometric standard star 


G191-B2B; and remove telluric features from the red chan- 
nels using the spectro-photometric standard. The data from 
the three channels were combined into a single spectrum and 
re-binned to a common scale of Ad = 0.85 A pixel !. An af- 
terglow spectral energy distribution (SED; $3) was then sub- 
tracted from the spectrum and the absolute flux calibration 
was fine-tuned using the late-time g’, 7', i’, z MODS and 
the Large Binocular Cameras (LBCs; Giallongo et al. 2008) 
photometry (Table 1), which is free from SN and afterglow 
emission (see §2.4). 


2.2. Optical and Near-Infrared Imaging 


The optical follow-up was performed with LBT in the 
Sloan filters g’ and r’ bands using the two twin MODS in- 
struments, and in the Sloan r’ filter with the Device Opti- 
mized for the Low Resolution (DOLoRes) optical imager and 
spectrograph mounted on the Telescopio Nazionale Galileo 
(TNG) on La Palma (Canary Islands, Spain). We have also 
observed GRB 200826A twice with AZT-22 1.5m telescope 
equipped with the SNUCAM CCD camera (Im et al. 2010) 
at the Maidanak Astronomical Observatory (MAO, Uzbek- 
istan). 

All MODS observations were taken under good seeing 
conditions (075 — 10). Due to the binocular capabilities of 
LBT, the observations in g' and r’ were obtained simulta- 
neously with both MODS instruments. During summer of 
2021, with the aim of expanding our host-galaxy photometry 
data set, we obtained further Ur’i’z’ and Rc imaging with the 
LBCs under variable seeing conditions (0/9 — 1” 5), followed 
by further MAO observations in the Rc band. LBT/MODS 
and LBC imaging data were reduced using the data reduction 
pipeline developed at INAF - Osservatorio Astronomico di 
Roma (Fontana et al. 2014) which includes bias subtraction 
and flat-fielding, bad pixel and cosmic ray masking, astro- 
metric calibration, and coaddition. We have also downloaded 
the Gemini/GMOS data at 28 and 46 d reported in Ahumada 
et al. (2021), and afterwards reduced them using the dedi- 
cated Gemini pipeline DRAGONS (Labrie et al. 2019). Fi- 
nally, we have downloaded the r’-band Gran Telescopio Ca- 
narias (GTC) data at 3.99 d reported in Zhang et al. (2021)*. 
TNG/DOLoRes, GTC and MAO images were reduced in a 
standard manner using PyRAF/IRAF tasks (Tody 1993). 

The first H-band observations were obtained using the 
LBT Utility Camera in the Infrared (LUCI, Seifert et al. 
2003) imager and spectrograph in conjunction with the Sin- 
gle conjugate adaptive Optics Upgrade for LBT (SOUL; 
Pinna et al. 2021), the 2nd generation adaptive optics (AO) 
system at LBT. Data were obtained 37.1 d after the GRB, 
close to the expected maximum light of the possible asso- 


3 PIL. Singer, program GN-2020B-DD-104. 
^ PI A. J. Castro-Tirado, program GTCMULTIPLE4B-20A 
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Figure 1. Left: 1' x 1’ field-of-view (FoV) finding chart derived from the TNG r’-band imaging obtained on 2020-09-27 UT. The star used for 
AO observations with SOUL and the LBT/LUCI camera is highlighted. The rectangular region is centered on the location of the afterglow, and 
its angular size is the same as the panels on the right. Right: Results of the H-band follow-up using the second generation SOUL Adaptive 
Optics with LUCI-1 at the Large Binocular Telescope. The top row shows the location of the host and the GRB (see the region in the left panel). 
For comparison, the bottom row shows the location of a star in the field. From left to right: A) the LUCI-1 AO H-band imaging obtained on UT 
02-10-2020 (1.e., close to the peak of the SN); B) LUCI-1 AO H-band image obtained four months later when only the host is visible; and C) 
the subtraction between these two epochs, where the transient is highlighted with a circle (top) and the field star is cleanly subtracted (bottom). 
The FoV of each panel is 276 x 2/0. The LUCI-1 AO FoV is 30" x 30". In each panel North is up, East is to the left. The transient lies at a 
projected distance of 0.75 kpc from the center of its host galaxy as discussed in $3.5. 


ciated SN. Compared to J and K bands, the H band was 
preferred for three reasons: 1) lower sky and thermal back- 
ground with respect to the K band; 2) the AO performance is 
better than in J, allowing the system to achieve a full-width at 
half-maximum (FWHM) < 0 1; 3) there are no strong emis- 
sion lines from non-stellar sources or from star-formation re- 
gions (like Ha in the J-band at the redshift of our target), and 
the host continuum is dominated by light from old, late-type 
stars. A relatively bright foreground star (R = 15.06 mag) 
11"4 away from the GRB location (and within the field-of- 
view) was used as the “AO Reference star" for the AO sys- 
tem. The AO system achieved an average FWHM of the 
point-spread function (PSF) ~ 0713. Subsequent observa- 
tions were able to achieve a FWHM ~ 0708 — 0/713 (see $2.3 
for more details). The observations were repeated on De- 
cember 1st (2020) and February 7th (2021), 97 and 158 d 
respectively, after the GRB event. 

In September 2021 we also obtained J-band imaging of 
the host with LBT/LUCI. All NIR data have been reduced 
using the same data reduction pipeline as for LBT MODS 
data ($2.2), which in case of LUCI includes dark subtrac- 
tion, flat-fielding, bad pixel and cosmic ray masking, and sky 
subtraction. Finally, the images were aligned to match the 
position of the AO reference star and co-added using IRAF 
tasks. The astrometry was calibrated against field stars in the 
GAIA DR2 catalogue (Gaia Collaboration et al. 2018) and 
has an astrometric precision of 0715. 

All data were analyzed by performing aperture photometry 
using DAOPHOT and APPHOT under PyRAF/IRAF. For the 


optical data, the photometric calibration was performed us- 
ing a set of field stars observed during the fifth MODS epoch 
(which had photometric conditions) and calibrated against 
the SDSS DRI2 catalog magnitudes of brighter nearby stars 
(Alam et al. 2015). NIR photometry was calibrated against 
the AO reference star and a pair of 2MASS stars observed 
immediately after the GRB field during the first epoch. The 
photometric calibration has a standard deviation of 0.02 
mags. To err on the side of caution, 2x this value was used 
as the photometric uncertainty. 


2.3. Image subtraction and photometry of the transient 


To remove the contribution of the host to the afterglow/SN 
photometry, a reference frame free of any transient light is re- 
quired. To better investigate any possible contamination from 
SN light in these earlier observations, the GRB 2008264 data 
sets were compared with GRB 980425/SN 1998bw (Galama 
et al. 1998), which is the standard template for GRB-SNe. 
Following Zeh et al. (2004) and Klose et al. (2019), the 
observed SN 1998bw light curves were first corrected for 
Galactic extinction, moved to rest-frame, interpolated with 
a polynomial curve at steps of 0.5 d and finally shifted to the 
GRB 2008264 observed bands, 7’, i’, and H, accounting for 
filter band differences and cosmological k-corrections. 

At the time of the last MODS r’ observation (78.9 d after 
the trigger) any SN 1998bw-like GRB-SN would have had 
r' ~ 27 mag, more than 2 mag fainter than its peak magni- 
tude (r' ~ 24.7 mag) at 20 — 40 d, when the first observations 
had been obtained. Therefore, it is acceptable to consider 


Table 1. Photometry of the Transient and Host Galaxy. 
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At? Magnitude Filter Telescope — Ref.^ 
(d) AB? 

0.2112 20.86+0.04 g ZTF [1] 
0.2304 20.70+0.05 r ZTF j 
0.2779 20.96+0.16 g ZTF d 
0.7663 > 20.2 Rc Kitab [2] 
1.4458 22.75+0.26 g ZTF [1] 
1.2125 > 21.3 r ZTF i 
1.2875 > 21.2 g ZTF d 
1.7351 > 20.4 Rc Kitab [2] 
2.1834 23.39+0.16 r LBT/MODS [3] 
2.1834 24.08+0.19 g' LBT/MODS ^" 
2.27 > 23.3 r FTN+LCO [1] 
2.28 > 23.4 g’ FTN+LCO i 
3.23 24.46-0.12 r Lowell ^ 
3.990 24514014 r GTC [3.4] 
7.21 > 23.5 r — LBT/MODS [3] 
7.21 > 23.6 g’  LBT/MODS ” 
15.7 > 24.3 Rc MAO » 
28.3 24.53+0.21 7 Gemini [1,3] 
28.3 > 25.6 r Gemini [1,3] 
31.9 > 25.3 r TNG [3] 
33.3 > 25.7 r  LBT/MODS ^" 
33.3 > 25.8 g’  LBT/MODS ” 
37.1 24.06 + 0.20 LBT/LUCI g 
46.1 25.36+0.26 7 Gemini [1,3] 
46.1 > 25.5 r Gemini [1,3] 
50.88 > 25.6 r  LBT/MODS [3] 
50.88 > 26.2 g’  LBT/MODS ” 
97.0 > 24.7 H | LBT/LUCI di 
74.8 22.85+0.11 r TNG » 
78.9 22.83 +0.04 r”  LBT/MODS ” 
78.9 23.12+ 0.04 g’ LBT/MODS ” 
158.5 22.84+0.20 H  LBT/LUCI i 
287 23.36+0.10 U*  LBT/LBC 7 
287 22.92+0.08 7 LBT/LBC ^ 
287 22.84 + 0.07 Rc LBT/LBC ae 
287 22.58+0.05 7 LBT/LBC » 
287 2225x011 2’ LBT/LBC r 
340 22.85 + 0.07 Rc MAO di 
382 23.29+0.08 B MAO » 
391 22.65+0.14 J LBT/LUCI d 


a Mid-time after the Fermi/GBM trigger in the observer frame. 


the late r’ MODS images free from any contamination by SN 
light and they can be used as templates for image subtraction. 
They can also be used to measure the flux of the host galaxy 
(^ 4 mag brighter in 7’ band) instead of the later LBC im- 
ages, which have a lower signal-to-noise ratio. These LBC 
images were instead used as reference for the image subtrac- 
tion of the 7’ imaging at ~ 50 d. Similarly, in the H band? 
any SN1998bw-like GRB-SN contribution at 160 d (H ~ 26 
mag at the time of our last LBT/LUCI observation) would 
have been ~ 2 mag fainter than its peak magnitude (H ~ 24.3 
mag). Differences of ~ 2 mag translate to a contribution from 
the SN to the late-time image of ~ 0.2 mag. This contribu- 
tion does not impact the analysis after the image-subtraction, 
but is included in the photometric errors. On the contrary, 
any SN1998bw-like GRB-SN would fade less than 1 mag be- 
tween 97 and 160 d, thus the last H —band observation (97 d) 
may also be affected by over-subtraction, and hence a photo- 
metric value for the transient that is too faint at 97 d, though 
the SN was likely too faint to contribute substantially at this 
point. For completeness we report the resulting upper limit. 

For the i’ band we note that at the time of the final Gem- 
ini observation (at 75 d, see Ahumada et al. 2020b) any 
SN1998bw-like GRB-SN would have had i’ ~ 25.5 mag, just 
~ 1.5—2 mag fainter than the peak magnitude at ~ 30 d (~ 24 
mag) which could be acceptable. Instead, this difference is 
certainly less in the case of the Gemini observation at 46 d 
(see Table 1 and Fig. 7). Therefore, we caution that using the 
value at 75 d as pure host can produce an over-subtraction, in 
particular at 46 d. Therefore, in the image subtraction of the 
i’band we use the LBC imaging obtained at 287 d. 

Before applying image subtraction, input and reference im- 
ages were aligned using the WCSREMAP package?. Im- 
age subtraction was then performed using a routine based on 
HOTPANTS’ (Becker 2015). This algorithm matches the PSF 
and count flux of both input images. The PSF is modeled via 
Gaussian functions in sub-regions of the original image. The 
software outputs a noise map of the resulting difference im- 
age, which is used to derive the uncertainties in the measured 
fluxes. The routine creates a grid of input Gaussian widths 
and sub-region sizes and searches for the solution with the 
lowest noise. 

In the case of the H-band imaging, the average FWHM 
during the first epoch was 0713 as measured from the AO ref- 
erence star in the LUCI-1 FOV? and was constant for nearly 
the entire observation. The small amount of data with larger 


P The photometry is not corrected for Galactic extinction. The first 
part (above the double line) is the result of image subtraction and 
represents the pure transient. The second part has been used to 
study the host and as references for the subtraction. 

© [1] Ahumada et al. (2021); [2] GCN 28306 (Belkin et al. 20202); 
[3] this work; [4] Zhang et al. (2021). 

d These are LBT/MODS acquisition images obtained using a ~ 2’ x 
2’ FOV instead of the full ~ 6’ x 6’. MODSI was equipped with 
the g' filter and MODS2 with the 7’ filter. 

* We have used the U spec for its better sensitivity compared to the 
U gessel filter. 


5 To derive the H-band flux of SN 1998bw at the redshift of GRB 200826A 
we had to extrapolate the rest-frame Jc band data of SN 1998bw (e.g., Zeh 
et al. 2004; Klose et al. 2019). 


6 http://tdc-www.harvard.edu/westools 
7 https://github.com/acbecker/hotpants 
8 The PSF is best modeled with a Moffat profile. 
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FWHMs were not used for analysis. A total of 1180 s of ex- 
posure time with FWHM ~ 0713 was used for analysis. For 
the other epochs the FWHMs ranged from 0708 to < 015. 
Therefore, only the exposures with a FWHM of the AO ref- 
erence star < 0713 (to match the first epoch) were selected. 
This yielded a combined image with a total integration time 
of 2400 s and 2700 s for the second and last epochs, respec- 
tively. In particular, the AO reference star in the last epoch 
has a FWHM ~ (’10, and this was used as the template for 
image subtraction. In contrast to the optical images, only two 
stars (the AO reference star and another star close to it) were 
useful to map the PSF and flux. Unfortunately, this impacts 
the optimal use of HOTPANTS to scale the template image 
to the flux of the science image. Instead, a better result was 
obtained using the PSF-MATCH task under IRAF and scal- 
ing the images to have the same flux for the two stars. The 
image subtraction of the two H-band images shows a tran- 
sient located at RA (J2000) = 00^27" 08:55, Dec. (J2000) 
= +34°01'3872 with an uncertainty of 0715 and calibrated 
against the MODS images. This is consistent with the posi- 
tion reported by Ahumada et al. (2020a). 

All image-subtracted data have been analyzed using aper- 
ture photometry as detailed in $2.2. Table | provides a 
summary of all the photometry of the transient and the host 
galaxy including values collected from the literature. 

Apparent magnitudes were corrected for Galactic extinc- 
tion using the Cardelli et al. (1989) interstellar extinction 
curve, a total-to-selective extinction of Ry = 3.1, and a red- 
dening along the line of sight of E(B — V) = 0.058 mag 
(Schlafly & Finkbeiner 2011). 


2.4. Modelling of the light curve and SED 


Although the g’ r’ light curves are quite sparse, they can be 
modeled with a single power-law with a = 0.94 0.05. How- 
ever, this model implies that the late afterglow is brighter than 
the deep r’ upper limits obtained at ~ 30 d after subtraction of 
the host-galaxy contribution. Therefore, the data were mod- 
eled with a smoothly broken power-law (Beuermann et al. 
1999) (Fig. 2). The last detections in the g' and r’ bands 
are at x 2.2 and « 4 d (observer-frame). Using SN 1998bw 
as a SN template (see below for more details), the SN it- 
self would have been g' = 29.8 mag and r' = 26.9 mag at 
2.2 d and 4 d, respectively. Therefore, it can safely be as- 
sumed to not contribute significantly to the optical transient 
and can be ignored in the afterglow fitting. This yielded de- 
cay slopes a; = 0.43 + 0.55, a» = 1.59 + 0.20, with a break 
time fy = 0.47 + 0.24 d, break smoothness n = 10 fixed. The 
post-break decay slope easily accommodates the magnitude 
upper limits at later times, with values beyond r’ ~ 27.8 mag 
and g’ ~ 28.3 mag at ~ 30 d. Therefore, any late emission 
detected via subtraction is not from the afterglow. 
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Figure 2. Light curves of the afterglow and SN component of GRB 
200826A in the g’r’i’H bands, after subtraction of the host-galaxy 
contribution. The data have been modeled with a smoothly broken 
power-law for the afterglow, and with a SN 1998bw template red- 
shifted to z = 0.7486 for the SN phase. In the g’ and r' bands, 
we forced the individual SN light curve components to go through 
the most significant deep upper limit in each band to obtain upper 
limits on the potential SN luminosity. Note the g’ band is offset by 
1 mag for reasons of clarity. See Tab. 1, $2.4 and 83.3 for more 
information. 


To investigate whether the optical data are affected by lo- 
cal dust extinction along the line of sight, a spectral energy 
distribution (SED) was modeled spanning the optical to the 
X-ray bands. All Swift XRT data available were used, which 
have a mean arrival time of 65840 seconds (0.762 d), along 
with the simultaneous g’ and 7’ band obtained from the bro- 
ken power-law fit (see above) which is free from host con- 
tamination. The redshift was fixed to match that of the host 
galaxy and the Galactic hydrogen column density was set to 
Ny = 6.02 x 10? cm? (Willingale et al. 2013). The optical- 
to-X-ray SED is well-fit (y2/d.o.f. = 33/47) by a power- 
law with spectral index 6 = 0.79 + 0.03 without need for 
additional dust extinction and gas absorption to improve the 
fit. The negligible dust extinction is well-supported by the 
optical-to-X-ray spectral index Box > 0.5 (Jakobsson et al. 
2004; Rossi et al. 2012). The optical-to-X-ray spectral slope 
was then used to extrapolate the late-time brightness in the i’ 
and H bands. For the first i’ detection at 28 d, the brightness 
of the projected afterglow was 27.6 mag, 3.2 mag fainter than 
the detection. The H-band afterglow is H ~ 27.4 mag at 37 
d, 3.4 mag fainter than the detection. 

To consider also the contribution of an emerging SN com- 
ponent, we have followed the work of Zeh et al. (2004) and 
Klose et al. (2019), and modeled the afterglow and SN com- 
ponent in a joint g’r’i’H fit (rest-frame far-UV and U Bz' 
bands). The analytical light curve of SN 1998bw and the 
k, s parametrization were used (see §2.4). The parameter k is 
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the luminosity factor and describes the observed luminosity 
ratio between the GRB-SN at peak time and the SN template 
in the considered band. The parameter s is a stretch factor 
with respect to the used template’. The observer-frame g’ 
band (rest-frame far-UV) is not covered by observations in 
the case of SN 1998bw (in contrast to the observer-frame r’ 
band, which is the rest-frame U band), so we assume that the 
flux density of the SN falls as v? to extrapolate and derive 
the corresponding far-UV SN 1998bw template (e.g., Klose 
et al. 2019). Note that the g’ and r’ data points in Figure 2 
after 2.2 and 4 d (observer-frame) are upper limits only. In 
order to constrain the flux density of a possible SN contri- 
bution in these bands, we considered these limits as proxy 
detections in the modelling. We also remind the reader that 
the i’- and H-band afterglow light curves have been spectrally 
extrapolated (see above) because there are no observations in 
these bands during the afterglow phase. The sparse data in 
i'H preclude determining the precise temporal evolution of 
the late emission, but do not preclude determining more in- 
formation on the luminosity. The stretch factor was fixed at 
s = | for all bands and the luminosity factor k was allowed 
to vary freely and individually in each band. The results of 
this joint fit are presented and discussed in $3.3. 


3. RESULTS 


A summary of the photometry of the transient resulting 
from the host-subtraction analyses is presented in Table 1. 
NIR emission (rest-frame ~ 0.93 um) is clearly detected 37 d 
(21.22 d rest frame) after the burst (see also Fig. 1). Thanks 
to our late reference images obtained with the LBCs, which 
are free from SN light, we also detect the bump in the i’-band 
at 28.3 and 46.1 d (14.9 and 26.3 d rest frame). Ahumada 
et al. (2021) report a bump that is ~ 1 magnitude fainter and 
they do not find the transient in the second epoch. This last 
detection allows us to better constrain the late evolution of the 
transient. In the r'-band, the presence of the bump could only 
be constrained beyond the 3c magnitude limits of r’ > 25.5 
mag, r' > 25.3 mag, and r’ > 25.7 mag in the first Gem- 
ini, TNG and LBT epochs, respectively (the results presented 
here supersede the preliminary results of Rossi et al. 2020a). 


3.1. GRB 200826A 


The duration of GRB 2008264, as reported by the differ- 
ent detectors and in different energy bands, is 1 — 2 s, placing 
it right at the divide between SGRBs and LGRBs in the dura- 
tion distribution commonly used to distinguish between these 
two classes (but see de Ugarte Postigo et al. 2011; Bromberg 
et al. 2012). 


? As in $2.3, the cosmological effects such as k-correction were considered 
when moving the SN template to the redshift of our GRB. For further de- 
tails on the method see Zeh et al. (2004), Klose et al. (2019), and references 
therein. 
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Figure 3. GRB 200826A (red) in the E,; — Eiso plane. GRBs with an 
associated SN are highlighted in green, outliers in blue (see §3.1). 
GRB data is from Amati et al. (2019). 
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Figure 4. The Too; — EH diagram for SGRBs (blue squares) 


and LGRBs (red circles) with corresponding cluster analysis results. 
6896 and 9596 confidence regions are shown by bold solid and thin 
dashed curves. Data from Minaev & Pozanenko (2020a,b). 


There are other events that lie in between the two popu- 
lations, but none has an observed T», shorter than 2 s and 
an associated SN which clearly points to a collapsar origin. 
For example, GRB 090426 (z = 2.609) has a duration shorter 
than 2 s but also has characteristics typical of LGRBs in- 
cluding softer spectra, a bright afterglow and a blue and low- 
mass host (Antonelli et al. 2009; Thóne et al. 2011; Nicuesa 
Guelbenzu et al. 2011, 2012, Kann et al. 2022, in prep.). It 
could be argued that one should consider the intrinsic du- 
ration in the rest frame (Ts9;). Doing so, the duration of 
GRB 200826A would be between 0.39 (Konus-Wind) and 
0.57 (GBM) s, i.e., well below 2 s. The only other GRB 
with an associated SN that could be classified as short be- 
cause of its intrinsic duration (Too; = 2.4 s, z = 0.86) is GRB 
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040924 (Wiersema et al. 2008; Soderberg et al. 2006; Huang 
et al. 2005). Many other cases exist at high redshift like GRB 
080913 (Greiner et al. 2009; Zhang et al. 2009), but the red- 
shift is too high to search for the associated SN (at least with 
the current facilities). 

Bromberg et al. (2012) showed that the 2 s duration com- 
monly used to separate collapsars and non-collapsars is in- 
consistent with the duration distributions of Swift and Fermi 
GRBs and only holds for old CGRO/BATSE GRBs. How- 
ever, the analysis of the Konus-Wind data (Svinkin et al. 
2020) showed that GRB 200826A has a duration typical 
of SGRBs. Thus, the short duration of GRB 200826A is 
instrument-independent. In this case Too alone is insufficient 
to determine which population GRB 2008264 is best asso- 
ciated with. Another feature typical of LGRBs is a non- 
negligible spectral lag (Norris et al. 1996). This parame- 
ter has been calculated and the details are provided in Ap- 
pendix A. We obtain a spectral lag of 96 + 38 ms, which 
lies in between the values found by Ahumada et al. (2021) 
and Zhang et al. (2021). In agreement with their analysis, 
we conclude that the measured spectral lag is more typical of 
LGRBs. 

It is only when both the spectral properties of GRB 
2008264 and its location in the hardness-duration plane are 
considered that it appears to be more similar to "long" GRBs. 
However, if this is indeed the shortest “long” GRB associ- 
ated with a SN, then additional evidence is needed to demon- 
strate that it does not belong to the SGRB class. Such a result 
would challenge the standard paradigm for SGRB and LGRB 
progenitors. 

In this respect, a very useful tool is the location of a GRB 
in the Ey; — Eiso plane, where LGRBs follow a strong cor- 
relation known as the “Amati Relation" (Amati et al. 2002; 
Amati 2006). In comparison, SGRBs populate a different 
region in this plane (e.g., D'Avanzo et al. 2014; Minaev & 
Pozanenko 2020a). Therefore, the fluence and spectral pa- 
rameters reported by the Konus-Wind and Fermi/GBM teams 
were used (Ridnaia et al. 2020; Svinkin et al. 2020; Mangan 
et al. 2020) to compute the values of Ej; and Eiso, and their 
uncertainties, for GRB 2008264. As can be seen in Fig. 3, 
this event is fully consistent with the Ep; — Eiso relation fol- 
lowed by LGRBs. This further supports the hypothesis that 
it is actually a "long" GRB with a duration in the short tail of 
the distribution. 

Recently, the energy-hardness (EH) parameter — a combi- 
nation of Ep,i and Eiso (EH = Eus — was introduced 
in Minaev & Pozanenko (2020a). Although GRB 200826A 
is placed inside the 95% cluster region of long bursts in the 
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Figure 5. The optical afterglow of GRB 200826A (red line) com- 
pared with the sample of extinction-corrected afterglows shifted to 
z = 1 from Kann et al. (2010, 2011, 2022, in prep.). Hereby, time 
and magnitudes are given in the observer frame, but assuming all 
GRBs are at z = 1 in a perfectly transparent universe. Light gray are 
LGRBSs, thicker black lines SGRBs with secure redshifts (solid lines 
connect detections, dashed lines connect upper limits). All magni- 
tudes are Vega magnitudes. The afterglow of GRB 200826A lies 
intermediately between those of secure LGRBs and secure SGRBs. 


Too; — EH diagram (Fig. 4), it is evident that this does not 
permit an unambiguous classification of a burst being close 
to the cluster region intersection (see also Zhang et al. 2021 
for a similar analysis and conclusion). Also the position in 
the Ej; — Eiso plane may not be sufficient for classification. 
Indeed, interesting examples such as GRB 201221D!! (Agiif 
Fernández et al. 2021) exist which follow the E; — Eiso rela- 
tion for LGRBs but are in the SGRB region in the Too; - EH 
diagram. 


3.2. The afterglow of GRB 200826A 


A well known-result is that most SGRB afterglows are 
much fainter than LGRB afterglows (Gehrels et al. 2008; Ny- 
sewander et al. 2009; Kann et al. 2010, 2011). Following 
the method of Kann et al. (2006), using a spectral slope of 
p = 0.79 and no line-of-sight extinction, in Fig. 5 the opti- 
cal afterglow of GRB 200826A was compared with a large 
sample of SGRB and LGRB optical afterglows (taken from 
Kann et al. 2006, 2010, 2011; Agüí Fernández et al. 2021, 
Kann et al. 2022, in prep.). Similarly, in Fig. 6 the X-ray 
emission of GRB 200826A was compared with that of ~ 500 


10 Epi = Ep,obs(1 + z) is the rest-frame photon energy at which the vF, spec- 
trum peaks, and Eiso is the isotropic-equivalent radiated energy as measured 
in a “bolometric” band, usually 1 keV — 10 MeV in the rest frame. 


!! With a rest-frame Fermi/GBM Too = 0.068 s (Hamburg et al. 2020), and 
Swift/BAT Too = 0.078 + 0.0196 s (Krimm et al. 2020). 
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Figure 6. The X-ray afterglow of GRB 2008264 (red line) in 
the context of the X-ray afterglows of 449 Swift GRBs with known 
redshifts. This event could be followed-up by Swift/XRT only after 
10ks, but the few data points obtained show that the luminosity of 
GRB 2008264 is in between the long and short populations. The 
grayscale on the right is used to convert a given luminosity and time 
into a fraction of bursts. 


GRB afterglows (of which 38 are SGRBs) detected in at least 
two epochs and with known redshift. The X-ray light curves 
where obtained from the Swift Burst analyzer'? (Evans et al. 
2010) and put in the same frame using the method described 
in Schulze et al. (2014). The afterglow of GRB 200826A 
is among the faintest afterglows of LGRBs but also among 
the brightest ones of SGRBs , both in the X-ray and optical 
bands. Thus, the afterglow luminosity is not a good method 
for discriminating between, nor understanding the nature of 
this burst. 


3.3. The late bump as supernova emission 


Although the sparse data precludes determining the evo- 
lution of the late bump seen in the transient light curve, a 
comparison can be made with other GRB-SNe to look for 
similarities. This comparison includes: GRB 980425/SN 
1998bw (Galama et al. 1998) GRB 060218/SN2006aj (Fer- 
rero et al. 2006); GRB 100316D/SN1010bh (Olivares E. et al. 
2012); GRB 130702A/SN 2013dx (D’Elia et al. 2015; Toy 
et al. 2016; Volnova et al. 2017; Mazzali et al. 2021); and 
GRB 161219B/SN 2016jca (Cano et al. 2017b; Ashall et al. 
2019b). In Fig. 7 we show only those GRB-SNe with com- 
parable peak luminosities and colors, including SN 1998bw 
which has been used to also fit extreme events such as SN 
2011kl (Greiner et al. 2015; Kann et al. 2019) and the SN 
bump of GRB 1405064 (Kann et al. 2021). Due to the non- 


12 https://www.swift.ac.uk/burst_analyzer/ 


negligible redshift of GRB 200826A, the observed r’-, i’- and 
H-band data sets correspond approximately to the rest-frame 
U, B, and z'-bands, respectively. In the following compari- 
son, we will refer to these rest-frame bands. 

The SN associated with GRB 2008264 appears fainter 
than SN 1998bw in the rest-frame B-band than SN 1998bw, 
but is of comparable luminosity or slightly brighter (x 0.5 
mag) than SN 1998bw in rest-frame z-band. The most strik- 
ing feature is the suppression in the rest-frame U-band, at 
least 1 mag compared to SN 1998bw at 28 d. This prompted 
a further investigation into GRB-SNe with redder colors, 
and produced a possible agreement with GRB 130702A/SN 
2013dx. This GRB-SN (like GRB 161219B/SN 2016jca and 
possibly GRB 050525 A/SN 2005nc, Della Valle et al. 2006) 
is characterized by a rise-to-maximum time At <~ 15 d in 
the respective rest frame (e.g., Belkin et al. 2020b), which is 
a parameter related to properties of the progenitor stars (e.g., 
González-Gaitán et al. 2015). 

In section 2.4 we have modelled the late light curve con- 
sidering the analytical light curve of SN 1998bw and the k, s 
parametrization following the work of Zeh et al. (2004) and 
Klose et al. (2019), which allows us to place quantitative con- 
straints on the above analysis. We remind the reader that 
there are no detections in g’r’ after 2.2 and 4 d (observer- 
frame), thus with our modelling we can only place deep up- 
per limits on the flux density of the SN contribution in these 
bands (see Fig. 2). In terms of luminosity, a potential SN 
contribution in the rest frame far-UV (observed g’ band) is 
found as ky < 1.5, an unremarkable limit to its luminosity. 
Therefore, although the spectral slope of the extrapolation 
mentioned above is insecure, this does not have any signif- 
icant impact on our modeling. However, there is clear UV 
suppression in the rest frame U-band, we find ky < 0.43, im- 
plying that the SN associated with GRB 200826A is less than 
half as luminous as SN 1998bw in the same rest-frame band. 
In the rest frame B-band the suppression is also strong com- 
pared to SN 1998bw with k; = 0.56 + 0.10. The rest-frame 
z’-band detection (observed H band) is, however, marginally 
brighter than SN 1998bw (ky = 1.18 + 0.22). Moreover, it 
lies close to the SN peak for s = 1, and thus a significantly 
faster or slower SN would be even brighter. When compared 
with the list of k, s values of Klose et al. (2019) and Cano 
et al. (20172) the k factor ranges between 0.5 and 2. This 
provides plausible support for a SN. 


3.4. Possible Kilonova emission 


In Fig. 8, the late-time photometry of the transient is com- 
pared with the light curve of AT2017gfo, the KN associ- 
ated with the binary NS merger GW 170817 and the SGRB 
1708174 (Abbott et al. 2017c). The KN light curves are ob- 
tained from observed spectra shifted and scaled to the same 
redshift as GRB 2008264, following the method described in 
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Figure 7. Light curve in i’ and H bands after removing the host component via image subtraction (see Table 1). Only data after 10 d (observer 
frame) are shown. The data are corrected for Galactic extinction (see § 2.3). The SN light curve templates obtained from SN 1998bw, 2013dx, 
and SN 2016jca are shown for comparison and to stress the large range of variability and color of GRB-SNe (See 83.3). Downward-pointing 


triangles are upper limits. 
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Figure 8. Light curves of GRB 200826A after host subtraction 
compared with the light curves of AT2017gfo (solid lines) from 
Rossi et al. (2020b) and with theoretical KN light curves boosted by 
a magnetar from Perego et al. (2017) with three different magnetic 
field strengths (dotted, dashed, dot-dashed lines; §3.4). Downward- 
pointing triangles are upper limits. 


Rossi et al. (2020b). The comparison shows that the luminos- 
ity of the transient was at all times at least 3 mag brighter than 
the peak luminosity of AT2017gfo. This is not extraordinary, 
because KNe z 10 times more luminous in the optical have 


been claimed (e.g., Gao et al. 2017; Fong et al. 2021). How- 
ever, what makes this event very different from all KNe ob- 
served so far is its slower temporal evolution: assuming that 
the rest-frame B-band detection (i’ band in observer-frame) 
detection is close to maximum light, i.e. ~ 16 d in the rest 
frame, then the peak would be more than ~ 15 d later than 
that of AT2017gfo. The analysis below investigates whether 
such behavior is theoretically possible. 

Modeling of KNe powered by radioactive decay of r- 
process nuclei showed that KNe produced by nearly equal- 
mass binary NS mergers are expected to have peak times and 
magnitudes similar (within a few d and one or two mag- 
nitudes, respectively) to the ones observed in AT2017gfo, 
especially in the case of face-on mergers characterized by 
lanthanide-free high-latitude ejecta (e.g., Radice et al. 2018b; 
Kawaguchi et al. 2020). Viewing angle and composition ef- 
fects possibly decrease the KN luminosity compared with 
face-on configurations (Korobkin et al. 2021). In the case 
of very unequal-mass binary NS or black-hole NS mergers, 
the peak evolution can be slower, especially in the rest-frame 
NIR bands, but with similar peak magnitudes (e.g., Barbieri 
et al. 2019; Bernuzzi et al. 2020; Kawaguchi et al. 2020; Zhu 
et al. 2020). These theoretical results seem to also disfavour 
the interpretation of the observed transient as a KN powered 
by radioactive decay. 

However, it is known from several studies (e.g., Gompertz 
et al. 2018; Ascenzi et al. 2019; Rossi et al. 2020b; Rastine- 
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jad et al. 2021), that the distribution of luminosities covers 
a large range, and some KNe can be ten times more lumi- 
nous than AT2017gfo, such as in the cases of those associ- 
ated with GRBs 060614 and 050724 (Gao et al. 2017). Even 
more extreme cases that are ~ 100 times more luminous have 
been claimed in association with GRB 070714B (Gao et al. 
2017) and more recently GRB 200522A (Fong et al. 2021; 
O’Connor et al. 2021), although the evidence is not very 
strong. To achieve this large blue luminosity, the KN can 
be boosted by energy deposition from a magnetar remnant. 
Indeed, long-lived merger remnants could have a significant 
impact on KN light curves. While larger amounts of ejecta, 
expelled due to angular momentum excess, can increase the 
peak magnitudes, times, and widths only marginally (e.g., 
Radice et al. 2018a), the presence of a fast-rotating magne- 
tar (in addition to being a possible GRB central engine) can 
boost the KN luminosity via spin-down luminosity (Yu et al. 
2013; Metzger & Piro 2014; Siegel & Ciolfi 2016). 

To explore this scenario, the KN model presented in Perego 
et al. (2017) was used and extended to consider spin-down 
energy injection, as well as sources at cosmological distances 
in which the fluxes, the (thermal) spectra, and the observer 
time are corrected for the cosmological redshift. The com- 
puted KN light curves are shown in Fig. 8. For simplicity, a 
spherically symmetric ejecta model is considered for which 
0.02 Mo of matter expands homologously at 0.2c. Given that 
the adopted model foresees the presence of a long-lived mag- 
netar, an effective gray opacity « = 1 cm? g^, typical of 
lanthanide-poor ejecta, was assumed. Po = 0.7 ms was fixed 
as a typical initial rotation period (Radice et al. 20182), while 
the strength of the magnetic field B and the energy deposi- 
tion efficiency for the spin-down luminosity ey were varied 
to obtain a plausible agreement with the observations. A rea- 
sonable match for B ~ 10!^ G and &, = 0.03 are shown 
in Fig. 8. For the considered bands, larger magnetic field 
strengths (B = 10P G) result in lower, more rapidly declin- 
ing peaks, and a plausible match with the data would require 
&n 2 1. Weaker field strengths (B = 10! G), on the other 
hand, produce more luminous, slower evolving peaks, which 
plausibly match with the data for e x 0.001. The ejecta prop- 
erties also influence the light curves, but less significantly. 

The calculations suggest that the observed transient could 
be compatible with the late emission of a spin-down pow- 
ered KN for plausible values of the magnetic field and ther- 
malization efficiency. Unfortunately, all of the magnetically- 
boosted KN models fail to explain the red color of the tran- 
sient, and in particular the strong UV suppression observed in 
the observer-frame r’ band, although modeling greatly sim- 
plifies the ejecta, the emission properties, and in particular 
its opacity. But most importantly, during the first 10 d the 
theoretical KN is comparable or brighter than the observed 
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Figure 9. Le PHARE modeling discussed in $3.5 of the SED ob- 
tained from the final LBT/LBC URci'z, MAO B, LBT/MODS g'r', 
LBT/LUCI J, and LBT/LUCI+SOUL H band imaging (see Table 
1). The blue points indicate where the photometry would be placed 
according to the best-fit model and without contribution from emis- 
sion lines. 


Table 2. Emission lines and their measured fluxes corrected for 
Galactic extinction. 


Lines flux 
[1077 erg/cm?/s] 


[OI] 43727 12.6 + 0.4 
[Nell] 43869 2.1 + 0.4 
Hel 43889 1.2 + 0.4 
[NeIII] 43968 1.2 + 0.4 
Hó 44102 1.8 + 0.4 
Hy 44340 2.8 + 0.4 
Hf 44861 6.4 + 0.6 
[OIII] 44959 8.0 x 0.3 
[ONI] 45007 26.5 + 0.9 


The FWHM of the lines is ~ 6.1 A. The mean redshift of the emis- 
sion lines is z = 0.748577 + 0.000065. 


data without accounting for the afterglow. Therefore, the 
magnetically-boosted KN models are excluded. 


3.5. The host galaxy 


As first reported in Rothberg et al. (2020a,b), the con- 
tinuum of the host galaxy was clearly detected between 
0.38 um — 0.98 um. A redshift of z = 0.748577 + 0.000065 
was determined from the simultaneous detection of the fol- 
lowing emission lines: [O I], [O IH] and [Ne IN] dou- 
blets, He I 43889, and Balmer Hó 44101, Hy 44340, and 
Hj 44861 (Fig. 10). Fluxes of the lines were measured using 
the slinefit code (Schreiber et al. 2018) which also allows 
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Figure 10. LBT/MODS observed spectrum of the host galaxy of GRB 200826A binned by 2.5 À for presentation purposes. The detected 


emission lines are marked. 


correction of the Balmer-line fluxes for the underlying stellar 
absorption. All the measured fluxes are reported in Table 2. 
The absorption lines of the Mg 112A2796,2803 doublet were 
also detected at a lower redshift (z = 0.7462 + 0.0002) which 
may possibly imply that they are due to an intervening Mg II 
absorber along the line of sight (Vergani et al. 2009; Chris- 
tensen et al. 2017). The rest-frame equivalent width for both 
lines in the doublet is W,esr = 2.4 + 0.6 A. 

Using the derived redshift as a fixed input, the optical/NIR 
SED of the host galaxy was modeled (Table 1) with the code 
Le PHARE (Amouts et al. 1999; IIbert et al. 2006). The best 
fit (y7/d.o.f. = 2.2/9) points to a galaxy dominated by a 
young population (0.04*005 Gyr) with global dust reddening 
E(B-V) ~ 0.2 mag using the Small Magellanic Cloud (SMC) 
curves from Pei (1992) and low stellar mass log M./Mo = 
8.6 + 0.2 (Fig. 9). The B-band dust-corrected absolute mag- 
nitude is Mg = —20.01 mag and the star-formation rate is 
SFR = 13.0+10:9 Mo/yr. With this large uncertainty, this SFR 
is consistent with what is measured from the emission lines 
(see below). The dust extinction is in agreement with the 
broad range of values found for LGRB hosts, also in com- 
parison with the SFR (Hunt et al. 2014; Japelj et al. 2016). It 
is different from the afterglow-derived negligible extinction, 
indicating that the GRB sight-line might be crossing a rela- 
tively dust-free region (e.g., by being placed at the front side 
of the host galaxy). Respective to Ahumada et al. (2021), we 
find a mass that is a factor ten lower. The computed lower 


mass is most likely because the NIR observations were ob- 
tained sufficiently late to be free of any contributions from 
the transient. 

Using the Balmer decrement, and assuming the theoreti- 
cal values for the H:y/H8 and Hó/H ratios, in the absence of 
dust, for case B recombination at a gas temperature of 10^ K 
(Osterbrock 1989), and a Milky Way (MW) extinction law? 
(Pei 1992) with Ry = 3.1, we derive an E(B-V) of ~ 0.16 
mag and Ay ~ 0.5 mag. Using extinction corrected H6 and 
[O rJ fluxes and the prescriptions provided by Kennicutt 
(1998) and Savaglio et al. (2009), the resulting star-formation 
rate (SFR) of the host galaxy is SFR ~ 4.0 Mo/yr. The SFR as 
determined from the line fluxes is within 2c of that derived 
from SED fitting. The high SFR obtained from two separate 
methods implies a specific star-formation rate (SSFR, SFR 
per unit stellar mass) of 1-3x1078 yr7!, among the highest 
measured so far in the LGRB host population (e.g., Savaglio 
et al. 2009; Hunt et al. 2014; Perley et al. 2013; Vergani 
et al. 2015; Japelj et al. 2016; Schulze et al. 2018; Modjaz 
et al. 2020). For comparison, this sSFR is about two orders 
of magnitude higher than for the Large Magellanic Cloud. 
Moreover, if the derived SFR were constant, this would im- 
ply that the time required for the host to build the measured 


13 Note that the difference between the MW and the SMC or LMC extinction 


laws is small in the wavelength range of Hô to Hf. 
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stellar mass would be about 0.03—0.1 Gyr, which is consistent 
with the upper limit obtained from the SED best-fit. Using 
several diagnostic tools in Lamareille (2010), Trouille et al. 
(2011), and Marocco et al. (2011) the presence of an AGN, 
which may mimic a high SFR, can be confidently excluded. 

The diagnostic ratios of [O I], [O II], and Hf emis- 
sion lines were used to derive the gas metallicity in the host 
galaxy. Many calibrators have been proposed in the past (for 
a review, see Maiolino & Mannucci 2019); the method pro- 
posed by Curti et al. (2017) was adopted which computes the 
best-fit metallicity by minimizing Y? in the space defined by 
the different possible diagnostics given the available emis- 
sion lines. This yielded 12 + log(O/H) = 8.31 + 0.02. This is 
in agreement with the value of 12 + log(O/H) ~ 8.3 obtained 
from using the [Ne III]43869/[O II] line ratio (e.g., Bian et al. 
2018), a method which has more uncertainties. The resulting 
sub-solar metallicity (Z/Z> = 0.4) is consistent with the value 
derived for the LGRB hosts at similar redshift (Japelj et al. 
2016). However, we note that the redshift evolution of the 
mass-metallicity relation in normal star-forming galaxies de- 
rived by Bellstedt et al. (2021) shows that galaxies at z = 0.62 
or 0.8 and stellar mass similar to the GRB 200826A host have 
a typical metallicity of about 12 + log(O/H) = 7.85 or 7.80, 
respectively, significantly below our measured value. Taken 
together, these results indicate that the host galaxy of GRB 
200826A is characterized by a typical metallicity but larger 
than usual SFR and sSFR for a LGRB host. 

For the adopted cosmological parameters and using the last 
deep H-band observation, the half-light radius of the host 
galaxy is Ra = 0716 or 1.2 kpc at the redshift of the GRB. 
This size is consistent with LGRB hosts of similar masses 
(Kelly et al. 2014). The high angular resolution from the 
LUCI-AO observations also makes it possible to measure 
the GRB-SN position within the host. The transient has an 
offset of R = 0/1 from the center of the host, which corre- 
sponds to 0.75 kpc. This corresponds to a normalized offset 
R/R, = 0.62 which is also consistent with the vast majority 
of LGRBs (e.g., Bloom et al. 2002; Blanchard et al. 2016). 


4. DISCUSSION 


Additional support for the collapsar scenario was recently 
reported by Rhodes et al. (2021). By analyzing the radio and 
X-ray light curves, they found evidence of a sharp rise peak- 
ing at 4 — 5 d after the trigger in the radio bands. Within 
the LGRB/collapsar scenario, they propose that the peak is 
produced by the synchrotron self-absorption frequency mov- 
ing through the radio band, resulting from the forward shock 
propagating into a wind medium. Is our modeling of the 
afterglow and optical-to-Xray SED in agreement with their 
claim? The spectral slope of 6 = 0.79 we have found in $2.4 
favors the emission in both optical and X-rays to be below the 
cooling frequency, i.e., 6 = (p — 1)/2, with p = 2.58 + 0.06 


following Sari et al. (1998), a typical value for synchrotron 
emission in GRB afterglows. If the afterglow is in the post- 
break phase when we measure the SED, the spectral and de- 
cay slopes are difficult to justify within the standard afterglow 
theory. Therefore, assuming that the afterglow is still in the 
pre-break phase (spherical expansion) the temporal decay in- 
dex for emission below the cooling break is « = 3(p — 1)/4 
and a = (3p—1)/4 for an ISM and wind environment, respec- 
tively. Using a = 1.59 + 0.20, we obtain p = 3.12 + 0.27 and 
p = 2.45 + 0.27, respectively (Granot & Sari 2002), and thus 
we find that propagation into a wind environment could well 
explain both spectral index and light-curve decay, in support 
of the analysis of Rhodes et al. (2021). While steep, a pre- 
break decay of a = 1.59 is not unheard of, the sample of Zeh 
et al. (2006) contains two out of 16 well-defined light curves 
with o; % 1.75 followed by clear breaks and much steeper 
post-break decay slopes ay ~ 3. In such a case, the early 
shallower decay of the afterglow of GRB 200826A could 
represent an optical plateau (e.g., Dainotti et al. 2020, Ron- 
chini et al., in prep.). In this scenario the late shallow decay 
Observed in the X-rays can be explained as energy injection 
from a new-born magnetar (see also Zhang et al. 2021), but 
also by the receding cooling frequency after the passage of 
the wind termination-shock radius (e.g., Schulze et al. 2011). 

All arguments presented in this work point to a massive- 
star origin consistent with a collapsar event despite the very 
short duration. Several theoretical routes have been presented 
to explain this. In one of the scenarios discussed by Zhang 
et al. (2021) and by Ahumada et al. (2021) to explain the 
short duration of the burst, if the progenitor was a massive 
star it is possible that the total duration of the engine was 
longer than the rest-frame duration of 1 s, but somehow the 
jet could not emit gamma-rays for a longer time. For example 
the jet could spend a longer time penetrating the envelope of 
the star. Alternatively, the jet might not be highly relativistic 
for longer than 1 s, so that no bright y-rays could be produced 
after this time like, e.g., in the case of a newborn magnetar 
that initially injects a baryon-loaded, neutrino-driven wind 
(see also Ghisellini et al. 2020). Zhang et al. (2021) also 
note that the energy injection from the newborn magnetar 
could explain the late shallow decay observed in the X-rays. 
Other alternatives exist like the mechanism known as the in- 
homogeneous jet model, proposed by Yamazaki et al. (2004), 
in which the jet consists of multiple sub-jets which we ob- 
serve as short GRBs. Therefore, a fraction of SGRBs may 
come from collapsars. Peng et al. (2021) have recently pro- 
posed that the burst could have originated from the collapse 
of a Thorne-Zytkow-like Object (TZIO), which consists of a 
central NS and envelope formed after a coalescence with a 
white dwarf. They found the collapse of such a TZIO can 
naturally explain the short duration of GRB 200826A and 
the interaction between the disk wind and the ejected mate- 
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rial can explain the “supernova bump". However, this model 
needs to be proven with the updated photometry presented 
here, and in particular with our result that the SN is similar 
to other GRB-SNe. We have also demonstrated the impor- 
tance of Adaptive Optics (AO) to studying GRBs. The NIR 
follow-up observations presented here are the first to detect 
and pinpoint the location of a SN associated with a GRB, 
track the SN evolution, as well as determine the properties 
of the host galaxy. The last published attempts were more 
than 11 years ago and served only to provide information on 
the host galaxy of GRB 060418 (Pollack et al. 2009), and in 
case of GRB 070610 to determine that a NIR transient was 
actually the counterpart of a galactic source (Castro-Tirado 
et al. 2008), making that GRB association uncertain (Kasli- 
wal et al. 2008). The only other approach similar to that out- 
lined in this paper was the case of GRB 071003 (Perley et al. 
2008). In that case, AO observations made possible to de- 
tect the afterglow which was localized between two bright 
stars, but those observations were too early and the redshift 
too high to detect the associated SN (z = 1.60435, and ~16 d, 
respectively). The results presented here demonstrate clearly 
that AO can provide a powerful tool for studying GRBs. 


5. SUMMARY AND CONCLUSIONS 


GRB 200826A is a temporally short GRB at z = 
0.748577 + 0.000065 with a rest-frame duration of ~ 0.5 s, 
below the threshold of 2 s commonly used to separate SGRBs 
and LGRBs. To better understand the nature of this event, 
we initiated a follow-up campaign spanning a period of 1 
year which involved the LBT telescope in Arizona (USA), 
the TNG telescope on La Palma (Canary Islands, Spain) and 
the Maidanak Astronomical Observatory (Uzbekistan). Tak- 
ing advantage of the adaptive optics capabilities of LBT, we 
were able to obtain deep H-band observations between 37 
and 159 d (observer frame), corresponding to the rest-frame 
z' band at 21 to 91 d after the burst trigger. Image subtraction 
shows a faint transient within its host galaxy. Moreover, im- 
age subtraction of archival Gemini i’-band images with late 
reference observations obtained with the LBT revealed an op- 
tical transient which is less affected by over-subtraction than 
had been when reference images obtained at an earlier time 
after the GRB had been used. Finally, we were able to put 
strong upper limits on the UV rest-frame luminosity thanks 
to our LBT and TNG 7'-band observations. Our results show 
that: 


e Despite its short duration, this event is consistent with 
the Epi — Eiso "Amati" relation followed by LGRBs. 
The spectral lag is also more typical of LGRBs. 


e It was followed by a relatively faint optical and X-ray 
afterglow with a luminosity that lies in between those 
of LGRB and SGRB afterglows. 


e The evolution and color of the late bump is in good 
agreement with other GRB-SNe, and especially with 
the fast rising GRB 130702A/SN 2013dx. GRB 
2008264 is one of the cosmologically most remote 
GRB-SNe detected to date, close to the sensitivity limit 
of the present generation of 8 to 10 m class optical tele- 
scopes. 


The possible alternative scenario of a genuine SGRB 
followed by a KN like AT2017gfo is not supported by 
the different evolution and luminosity of the light curve 
of the observed transient. 


The host galaxy of GRB 2008264 is remarkable be- 
cause it is typical of an LGRB host galaxy, but with 
higher SFR and sSFR rates than expected. 


The GRB lies at a projected distance of 0.75 kpc from 
the center of its host galaxy, which is consistent with 
the majority of LGRBs. 


In summary, with the obtained results we are faced with 
a GRB of short duration which exploded in a star-forming 
galaxy, with a moderately faint afterglow, emitted by a jet 
most likely propagating into a wind environment, and fol- 
lowed by a bump in the light curve whose color and lumi- 
nosity are typical for a GRB-SN. Thus we firmly classify this 
burst as a collapsar event. This evidence, together with the 
analysis of the energetics of this burst, further weakens the 
effectiveness of simple duration as an indicator of the source 
of a GRB. In addition, strong support is provided to theoret- 
ical predictions that collapsar-produced events may have an 
observable duration well short of the classical short/long di- 
vide (about 2 s), and down to 0.5 s or less (e.g., Bromberg 
et al. 2013). 

In the next years future missions, like the Space Variable 
Objects Monitor (SVOM; Paul et al. 2011) and the Gamow 
Explorer (White et al. 2021) will offer a combination of ex- 
tended sensitivity and energy bands that will increase the 
number of known GRBs. This, when coupled with the im- 
proved capabilities of the new generation of extremely large 
telescopes will allow us to observe both the GRB and SN 
components resulting from collapsar explosions in increasing 
numbers and at higher redshifts (e.g., Maiorano et al. 2018; 
Rossi et al. 2018). In this scenario, the first detection of a 
SN with AO observations represents what ground-based tele- 
scopes can achieve: they will not just offer a sharper view of 
the GRB-SN location within its host, but have the necessary 
depth to discover GRB-SNe at larger redshift than what has 
been possible up to now from the ground, and yet at similar 
wavelengths in rest-frame. Eventually, these future facilities 
will assess whether peculiar events like GRB 2008264 are 
actually the result of the rich variety of the collapsar phe- 
nomenon (Amati 2021). 
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Facilities: Swift(XRT), LBT, TNG, MAO 


Software: IRAF (Tody 1993), PyRAF (Science Software 
Branch at STScI 2012), DRAGONS (Labrie et al. 2019), 
HOTPANTS (Becker 2015), Le PHARE (Arnouts et al. 1999; 
Ilbert et al. 2006), WCSTools (Mink 2019) 
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APPENDIX 


A. SPECTRAL LAG 


We first extracted the Fermi/GBM light curves with | ms resolution for the best exposed Nal detectors 6, 7, 8, 9, and 11 and 
calculated the corresponding total light curve in the same energy bands as Zhang et al. (2021): 10-20 keV, and 250-300 keV, 
respectively. We linearly interpolated the background counts of the 16-ms binned curves by fitting the time windows [—-10 : 
-2]U [7 : 10] s. 

For the spectral lag calculation, we then restricted the following analysis to the time interval [-0.5 : 4] s and calculated the 
cross-correlation function (CCF) between the background-subtracted profiles of the two energy bands, both using the original 
binning of 1 ms as well as that of 16 ms. A positive lag corresponds to the softer band lagging behind the harder one. The 
peak of the CCF was found to be about 100 ms. To provide a more accurate estimate, we fitted the CCF around the peak with 
a third-degree polynomial between —0.4 and 0.6 s and it gave 95 ms. To estimate the uncertainty, we carried out the following 
simulations: starting from the 16-ms profiles, we obtained smooth versions of these curves by adopting the L1 trend filtering by 
Politsch et al. (2020) properly adapted to model GRB light curves (e.g., see Ghirlanda et al. 2021). We then obtained 1000 random 
realisations of both profiles by assuming uncorrelated Poisson noise, assuming the total expected counts for each temporal bin 
(i.e., interpolated background plus smoothed GRB signal). For each pair of simulated profiles the corresponding CCF was fitted 
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Figure 11. Distribution of the spectral lag analysis. 


around the peak following the same procedure as for the real CCF. We finally collected the resulting distribution of 1000 lags 
shown in Figure 11: this is approximately normal with mean and standard deviation (96 + 38) ms, which therefore represents our 
estimate with uncertainty for the lag between the chosen energy bands. Our value looks slightly shorter than the analogous result 
by Zhang et al. (2021), who obtained (157 + 51) ms, even though statistically compatible and well in accord with values obtained 
by those authors for adjacent energy bands (see their Extended Data Figure 4). 

We finally repeated the same simulations using the 1-ms light curves, still adopting the same smoothed versions properly 


calibrated to 1-ms bins, and obtained a very similar result. 
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